Ising model with periodic pinning of mobile defects 
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A two-dimensional Ising model with short-range interactions and mobile defects describing 
the formation and thermal destruction of defect stripes is studied. In particular, the effect of a 
local pinning of the defects at the sites of straight equidistant lines is analysed using Monte Carlo 
simulations and the transfer matrix method. The pinning leads to a long-range ordered magnetic 
phase at low temperatures. The dependence of the phase transition temperature, at which the 
defect stripes are destabilized, on the pinning strength is determined. The transition seems to be 
of first order, with and without pinning. 
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I. INTRODUCTION 

Striped magnetic structures in high-temperature su- 
perconductors and related materials have attracted much 
interest for more than a decade, both theoretically and 
experimentally Q, Q, El 0, IB, Q • In that context, moti- 
vated by recent experiments on (Sr, Ca, La)i4Cu2404i 
0, S| a class of rather simple two-dimensional Ising 
models has been introduced describingthe formation and 
thermal destruction of defect stripes Q ■ 

The model consists of spin- 1/2 Ising variables, mimick- 
ing Cu 2+ ions, and non-magnetic defects, .5 = 0, corre- 
sponding to holes. The spins are arranged in chains with 
antiferromagnetic interactions, J a < 0, between neigh- 
boring spins in adjacent chains. Along the chains, neigh- 
boring spins are coupled ferromagnetically, J > 0, while 
next-nearest neighbor spins separated by a defect inter- 
act antiferromagnetically, Jo < 0. The defects are al- 
lowed to move along the chain through the crystal. The 
mobility of the defects is determined by the changes in the 
magnetic energy encountered during their motion (an- 
nealed Ising model). 

In a 'minimal variant' of the model, the couplings in 
the chains, J and |Jo|, are assumed to be indefinitely 
strong. The minimal model has been shown to describe 
the formation of defect stripes, oriented perpendicular to 
the chains, whose coherency gets destroyed at a phase 
transition. At the transition, one observes a pairing ef- 
fect for the defects in the chains, reflecting an effectively 
attractive interaction between defects mediated by the 
magnetic interaction between the chains, J a . The ther- 
mal behavior of the full model, choosing experimentally 
realistic values of the couplings in the chains, resembles 
closely that of the minimal model . 

The aim of this paper is to study the impact of a local 
defect pinning energy of strength E p on thermal proper- 
ties of the minimal model. In the experimentally stud- 
ied (Sr, Ca, La)i4Cu2404i compounds 0, 13], holes are 
pinned by Ca- or La-ions, which, in turn, are rather 
immobile. In the following, we assume that the fixed 
pinning sites form straight equidistant lines perpendicu- 
lar to the chains, with the number of pinning sites being 
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FIG. 1: Sketch of the interactions in the Ising model on a 
square lattice with periodic pinning of mobile defects. 



equal to the number of defects. Beyond the specific ex- 
perimental motivation, the model is hoped and believed 
to be of genuine theoretical interest. 

Of course, the model still allows for thermal fluctu- 
ations of the defect stripes at finite pinning strength. 
Indeed, the instability of the defect stripes and the ef- 
fects of the pinning on the spin ordering are intriguing 
features of the present model. In particular, at low tem- 
peratures spin correlations are expected to become long- 
ranged for non- vanishing pinning, while they decay alge- 
braically when E p = • The dependence of the phase 
transition, at which the defect stripes get destroyed, on 
E p is an interesting aspect of the model as well. Without 
pinning, the transition temperature had been estimated, 
but the type of the transition had not been studied. 

The layout of the paper is as follows. In the next 
section, we shall introduce the model and the methods, 
Monte Carlo simulations and transfer matrix calcula- 
tions. Results will then be presented and discussed in 
Sect. III. Finally, a short summary concludes the article. 



II. MODEL AND METHODS 

We consider an Ising model on a square lattice, set- 
ting the lattice constant equal to one. Each lattice site 
is occupied either by a spin, Sij — ±1, or by a 
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FIG. 2: Typical Monte Carlo equilibrium configurations of the minimal model, 6 = 0.1 and q p 
temperatures kBT/\J a \ = 0.8 (a), 2.3 (b), and 2.9 (c). Only parts of the systems are shown. 



1.0, of size L = M = 40 at 



defect corresponding to spin zero, Sij = 0, see Fig. 1. 
The defects are mobile along one of the axes of the lat- 
tice, the chain direction. The sites in the j-th chain are 
denoted by We assume a ferromagnetic coupling, 

J > 0, between neighboring spins, Sij and Si±i,j, along 
the chain, augmented by an antiferromagnetic interac- 
tion, J < 0, between those next-nearest spins in the 
same chain, which are separated by a defect. Spins in 
adjacent chains, Sij and Sij±i, are coupled antiferro- 
magnetically, J a < 0. Usually a minimal distance of two 
lattice spacings between neighboring defects in a chain 
is assumed, i.e. two defects are separated by at least 
one spin due to strong short range repulsion between de- 
fects (alternatively, one may introduce an additional fer- 
romagnetic coupling between spins separated by a pair of 
nearest-neighboring defects). A local pinning potential 
acts on the defects, lowering the energy of the defects at 
fixed sites by an amount E p . In the following, we choose 
pinning sites along equidistant straight lines, i = i v , per- 
pendicular to the chains with the number of pinning sites 
being equal to the number of defects, Nd. Accordingly, 
the Hamiltonian of the model may be written as 

Ti. = — *^^[JSijSi±i j + JoSijSi±2j{l — Sf±i > j) 

+J a SijSi y j±i + E p (l — Sfj)5 iiip ], (1) 

see Fig. 1. We assume that the number of defects is the 
same in each chain, determined by the defect concentra- 
tion O, denoting the total number of defects divided by 
the total number of sites, N#/N. In this study, we set 
O = 0.1, where the distance between the pinning lines is 
then ten lattice spacings. 

In the following we consider the 'minimal' variant of 
the model by assuming the couplings in the chain, J and 



| Jo I, to be indefinitely strong [jj. Thence the spins form 
intact clusters in the chains between two consecutive de- 
fects, and neighboring spin clusters have opposite sign. 
Thermal quantities depend only on, say, kBT/\J a \ and 
the ratio q p — E p /\J a \. 

To study the minimal model with pinning of mobile 
defects, we used Monte Carlo techniques and the 
transfer matrix method [TlT ]. 

In the simulations, a new configuration of spins and 
defects may be generated by exchanging a defect with a 
neighboring spin in a chain, reversing the sign of the spin 
to keep intact spin clusters. The energy change associ- 
ated with this elementary process is determined by J a 
and E p , see the Hamiltonian (1). As usual, the related 
Boltzmann factor determines the probability of accept- 
ing the new configuration Of course, simulations are 
performed on finite lattices with N = LxM sites, L being 
the number of sites in a chain. We shall present results for 
L = M. We employ full periodic boundary conditions. 
To investigate finite size effects, the linear dimensions, L 
and M, were varied from 20 to 320. Typically, runs of 
at least a few 10 6 Monte Carlo steps per defect were per- 
formed, averaging then over such realizations to estimate 
error bars. The pinning strength, q p = E p /\J a \, ranged 
from to 2.0. 

The transfer matrix calculations were done in the stan- 
dard way [ijj with the matrices representing the inter- 
actions of the entire chains. All eigenvalues and eigen- 
vectors were computed numerically, enabling us to derive 
quantities for arbitrary M, being finite or infinite. Study- 
ing the case O = 0.1, L was chosen to be 20, with two 
defects per chain. Larger systems, i.e. with L being at 
least 40, are outside the current reach of computer facili- 
ties. Of course, one may study the case of more than two 
defect stripes in the case of L — 20 by enlargening the 
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defect concentration. We shall consider here, however, 
only the case of a fixed value of 8 = 0.1. q p ranged from 
to 5.0. 

Physical quantities of interest include the specific heat, 
C, and spin correlation functions (depending, in general, 
on the distance from the pinning lines, i.e. on i), parallel 
to the chains, 



Gi(i,r)= \ Y J {Si, j S i+r , j )\ I L. (2) 
and perpendicular to the chains, 



G 2 (i,r)=\y^(S i , j S itj+r )j/L, (3) 

considering systems with M = L. Without pinning, the 
defects are expected to be delocalized so that there is 
full translational invariance, and the spin correlations do 
not depend on i. Note that in the thermodynamic limit 
for infinitely large distance, r — > oo, the perpendicular 
correlations G2(r) determine the profile of the squared 
magnetization 

m 2 (i)= lim m\{i) = lim G 2 (i,L/2) (4) 

L — >oo L — >oo 

We also calculated less common microscopic quanti- 
ties which describe the stability of the defect stripes and 
the ordering of the defects in the chains. In particular, 
we computed the average minimal distance, d m , between 
each defect in chain j, at position (idij), and those in 
the next chain, at (i' d , j + 1), i.e. 

d m = y^(min \i d - i' d \)/N d , (5) 

id 

dividing the sum by the number of defects, N d . Further- 
more, we calculated the cluster distribution, n d (l), de- 
noting the probability that consecutive defects in a chain 
are separated by I spins, in analogy to the distribution 
of cluster lengths in percolation theory Our main 

emphasis will be on pairs of defects with 1 = 1. Finally, it 
turned out to be quite useful to visualize the microscopic 
spin and defect configurations as encountered during the 
simulation. 



III. RESULTS 

In the ground state, T — 0, of the minimal model, the 
defects form straight stripes perpendicular to the chains, 
separating antiferromagnetic domains of spins. Without 
pinning, E p = 0, the ground state is highly degenerate. 
Each arrangement of defect stripes separated by at least 
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FIG. 3: Profiles of the absolute magnetization Mi(fe), at pin- 
ning strength q p = 0.2 (full symbols) and 2.0 (open symbols), 
with k — 1 (squares), 2 (diamonds), 3 (triangles up), 4 (tri- 
angles left), 5 (triangles down) and 6 (triangles right). Re- 
sults have been obtained from simulations of systems of size 
L — M = 160. 



two lattice spacings has the same lowest possible energy, 
resulting in an exponential decay of the correlations Gi 
parallel to the chains, while the spins are perfectly cor- 
related perpendicular to the chains @. By introducing 
the pinning potential, E p > 0, the defect stripes coincide 
with the pinning lines, at i = i p . Obviously, G\ contin- 
ues to oscillate, but now with a constant amplitude. Of 
course, the spin correlations perpendicular to the chains, 
G2(r), are equal to 1 for even distances r and —1 for odd 
distances r, when staying away from the pinning lines, 

Increasing the temperature, T > 0, the defects are al- 
lowed to move so that the stripes start to meander and 
finally break up, as exemplified in typical Monte Carlo 
configurations depicted in Fig. 2. Due to the pinning the 
defects tend to stick to the pinning lines at low tempera- 
tures. The detachment or depinning of the defects from 
those lines is expected to occur without phase transi- 
tion, as had been shown in the framework of SOS models 
with pinning The mapping of the minimal model 

onto the standard SOS model has been discussed before 
0. However, once the defects take positions far from 
the pinning sites, the magnetic interactions may mediate 
effectively attractive couplings between the defects. As 
for vanishing pinning 9], these couplings, absent in stan- 
dard SOS models, may eventually destroy the coherency 
of the defect stripes through a phase transition, as will 
be discussed below. We shall provide numerical evidence 
that the transition is of first order. The effect of the pin- 
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FIG. 4: Specific heat, C, at q p = (squares), 0.2 (circles), 0.5 (diamonds), and 1.0 (triangles), for systems of size (a) L = M = 
20, showing results from transfer matrix calculations (solid line) and simulations, and of size (b) L = M = 80, obtained from 
simulations. 



ning on the meandering and breaking up of the stripes, 
for various physical quantities, is exemplified in Figs. 3 
to 7. Note that in most of the figures we did not include 
error bars since they were, typically, not larger than the 
size of the symbols. Such a statement would not hold 
for appreciably shorter Monte Carlo runs because of the 
rather slow fluctuations of the defect stripes. 

At T > 0, without pinning, E p = 0, the model shows 
no magnetic long-range order. The spin correlation func- 
tion parallel to the chains, G\, has been shown, doing 
a free-fermion calculation, to decay algebraically at low 
temperatures 0. Indeed, our new Monte Carlo results 
both for G\ and G2 are consistent with such an alge- 
braic decay in the low-temperature phase characterized 
by meandering defect stripes whose positions can fluctu- 
ate rather freely. In particular, for finite systems of size 
L x L, the profile of the absolute value of the magnetiza- 
tion, |mi(i)| = yj m? L (i) , reflects the translational invari- 
ance, i.e. it does not depend on i, and it decreases signifi- 
cantly with increasing system size L. In marked contrast, 
with pinning, E p > 0, at low temperatures long-range 
magnetic order sets in, as seen easily from the profiles 
of the absolute magnetization between two pinning lines. 
The profiles are denoted in the following by Ml (A;) with 
k running from 1 to 11; k = 1 and k = 11 denote the two 
pinning lines, the center line in between them is at k = 6. 
Obviously, one has Ml (12 — k) = Mz,(fc) for reasons of 
symmetry. Examples of pertinent profiles are displayed 
in Fig. 3 at weak, q p = 0.2, and strong, q p = 2.0, pin- 
ning. Long-range order at low temperatures follows from 
the fact that the magnetization especially near the cen- 
ter between the two pinning lines is largely independent 
of system size. At high temperatures, the magnetiza- 



tion decreases appreciably with increasing system size, 
tending to zero in the thermodynamic limit. Indeed, 
finite-size analyses allow one to locate the phase tran- 
sition temperature as a function of the pinning strength, 
T c {% = E p /\J a \). Estimates agree with those obtained 
from analyses of the specific heat C, to be discussed next. 
Note that Mj,(fc) (or an average over these absolute line 
magnetizations) may be considered as the order param- 
eter of the problem. 

Results for the specific heat C are depicted in Figs. 
4a and 4b for lattices with linear dimension L = 20 and 
80 at pinning < q p < 1.0. At fixed pinning and vary- 
ing temperature, one observes two maxima in C . The 
maximum at the lower temperature is almost indepen- 
dent of the system size, and it stems from the mean- 
dering of the defects stripes with few excitations, i.e. a 
small kink density, as we checked by analysing and sim- 
ulating corresponding SOS or TSK (terrace-step-kink) 
models 0, Q, Ej| with pinning, similarly to the case 
without pinning !)|. The lower maximum is shifted to- 
wards higher temperatures when increasing the pinning 
strength E p . It may eventually be masked by the up- 
per maximum. The upper maximum of C, occurring 
at T max (L), signals the instability of the defect stripes 
due to thermally excited large fluctuations of the defect 
positions. At strong pinning, these fluctuations are ex- 
pected to set in once the defects start to detach in sig- 
nificant numbers from the pinning lines, giving then rise 
to a large specific heat, see Fig. 4b. In any event, the 
height of the second maximum increases clearly with in- 
creasing system size, indicating a phase transition in the 
thermodynamic limit, L — > 00. To estimate the transi- 
tion temperature, we plotted T max (L) versus X/L, with 
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L going up to 160, see Fig. 5. From a linear extrap- 
olation one may approximate the phase transition tem- 
perature T c (q p ) = T max (L = oo). T c (q p ) is found to 
increase monotonically with q p . More specifically, we ob- 
tain the following estimates from the data depicted in 
Fig. 5.: kBT c (q p )/\J a \ = 1.1 ±0.1 at q p = 0.2 (being close 
to the estimate at q p — [<j]), 1.30 ± 0.1 at q p = 0.5, 
1.55 ± 0.1 at q p = 1.0, and 2.10 ± 0.05 at q p = 2.0, with 
error bars reflecting some of the uncertainty in the linear 
extrapolation. Finite size analyses for other quantities 
lead to similar estimates for the possible transition tem- 
perature, as already mentioned in context of the magne- 
tization profiles. 

With pinning, the magnetization changes more and 
more drastically for larger systems close to T max (L), com- 
pare to Fig. 3. This behavior may suggest that in the 
thermodynamic limit the phase transition is of first or- 
der, with a jump in the magnetization at T c . To clarify 
this aspect, we determined the perpendicular correlation 
length, following from G2, when approaching T c from 
high temperatures. The correlation length may be esti- 
mated from analyzing the function [16] 

Ceff(r) = " {-n®-) 

with G 2 (r) = 5^ I G 2 (i,r) I /L (6) 

i 

Typically, the 'effective correlation length' £ e ff( r ) m ~ 
creases rather quickly monotonically for small r until 
it acquires a plateau-like behavior, and finally it rises 
steeply due to the finite size effect and periodic bound- 
ary conditions. Obviously, at a plateau of height £o, one 
has G2 oc exp(— r/£o)- Indeed, in the thermodynamic 
limit for T > T c , the height of the plateau at large r 
obviously corresponds to the standard correlation length 
£. Much care is needed close to the transition because 
very large system sizes may have to be studied to get 
an extended plateau. From simulations of systems with 
L = M = 160, we determined the correlation length ver- 
sus temperature, at various fixed q p . Using linear extrap- 
olation near T c (q p ), see above, we estimate the perpen- 
dicular correlation length at the transition. It is found to 
increase from about 20 lattice spacings at E p /\J a \ = 2.0 
to about 30 lattice spacings at E p = 0; i.e., it is finite. 
This finding supports the suggestion that the destruction 
of the defect stripes occurs through a phase transition of 
first order, with and also without pinning. A remark of 
caution may be added for the case of vanishing pinning. 
There, spin correlations in the low-temperature phase 
decay algebraically, and one might expect a transition 
of Kosterlitz-Thouless type. As has been noted before, 
however, algebraic order can be also destroyed by a tran- 
sition of first order fiH . 

The destruction of the defect stripes can be seen rather 
directly in the average minimal distance between defects 
in adjacent chains, d m . In Fig. 6, simulational data for 
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FIG. 5: Size dependence of the location of the maximum in 
the specific heat, T m&x (L), as obtained from simulations, at 
q v = 0.2 (squares), 0.5 (diamonds), 1.0 (triangles up), and 2.0 
(triangles left) for L = M ranging from 20 to 160. 



system sizes L — M ranging from 20 to 160, at q p = 2.0, 
are displayed. The temperature dependence of d m re- 
sembles closely the one found for the model without pin- 
ning (31 . While at low temperatures d m (T) does not de- 
pend significantly on the system size, it starts to rise 
rapidly at some characteristic temperature, correspond- 
ing to T max (L) in the case of the specific heat, with the 
height of the maximum in the temperature derivative of 
d m increasing strongly with larger system size. The lo- 
cation of the maximum, signalling the breaking up of the 
stripes, moves to lower temperatures as L gets larger. 
The quantitative behaviour is quite similar to the one of 
the specific heat and the magnetization profiles, for the 
various pinning strengths q p = 0, 0.2, 0.5, 1.0, and 2.0. 

The destabilization of the stripes seems to be driven 
by effectively attractive couplings between consecutive 
defects in a chain, mediated by the spin interactions J a 
(possibly reminiscent of the spin-bag mechanism [ijj]). 
Indeed, effectively attractive couplings may occur when 
two such defects, say, in chain j, at sites and 
(i + m,j), are displaced strongly with respect to cor- 
responding defects in adjacent chains, j ± 1, so that the 
spins in those chains at sites in between ± 1) and 
(i + m,j ± 1) have the same sign as the spins between 
the two defects in chain j. Such a situation may be re- 
alized, for instance, when three defects in chain j are in 
a cage of four defects in total, at, say, sites (i, j ± 1) and 
(i + k,j ± 1), in the neighboring chains with spin clus- 
ters of the same sign between the two pairs of defects 
in these chains j — 1 and j + 1. Then two of the three 
defects in the cage will move towards each other 0. In 
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any event, due to the effectively attractive coupling, me- 
diated by J a , two consecutive defects in chain j tend to 
form a pair of next-nearest neighboring defects having 
the minimal separation distance of two lattice spacings. 
The temperature and size dependence of the probability 
to find such pairs of defects, given by the pair probability 
rid(l — 1), is depicted in Fig. 7, at fixed pinning strength, 
q p = 2.0, and various system sizes. In general, the pro- 
nounced increase of the pair probability occurs close to 
the temperature T max (L), where other quantities signal 
the thermal instability of the defect stripes as well. For 
larger system sizes the increase in 71^(1) gets sharper and 
sharper in accordance with a transition of first order. At 
strong pinning, the pair probability rises quite drastically 
already in systems of moderate size, see Fig. 7, possibly 
reflecting the moderate correlation length at the transi- 
tion, as discussed above. 

Note that the type of stripe instability we observe here 
is not included in standard descriptions of wall instabil- 
ities in two dimensions jlTl 12(1 l2lL l22j. where either the 
number of walls is not fixed, giving rise to incommen- 
surate structures, or dislocations play an important role 
in the context of melting of crystals. Also the bunching 
of steps in TSK models with attractive step-step inter- 
actions or instabilities in polymer filaments due to 
attractive couplings 0, H3| are rather different from the 
destruction of defect stripes due to the pairing of defects 
induced by the inter-chain magnetic interactions J a . 



IV. SUMMARY 

In this paper a two-dimensional Ising model with peri- 
odic local pinning of mobile defects has been studied. Al- 
beit the model has been motivated by recent experiments 
on cuprates with low-dimensional magnetic interactions, 
the model is believed to be of genuine theoretical interest 
as well. 

In particular, based on Monte Carlo simulations and 
transfer matrix calculations, the model is found to de- 
scribe the pinning, meandering and, finally at higher tem- 
peratures, the destruction of defect stripes. 

The pinning gives rise to a long-range ordered mag- 
netic phase at low temperatures while magnetic corre- 
lations decay algebraically at low temperatures without 
pinning. 

The thermal instability of the defect stripes, which 
had been already identified for vanishing pinning, shifts 
towards higher temperatures as the pinning strength 



increases. The instability is signalled by pronounced 
anomalies, among others, in the specific heat, in the mag- 
netization profile, in the probability of defect pairs with 
shortest separation distance, and in the average mini- 
mal distance between defects in neighboring chains. The 
breaking up of the stripes is caused by an effectively at- 
tractive coupling between the defects mediated by the 
inter-chain interactions between spins in adjacent chains. 
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FIG. 6: Average minimal distance between defects in adjacent 
rows d m (T), at q p = 2.0, simulating systems of size L = M = 
20 (squares), 40 (diamonds), 80 (triangles), and 160 (circles). 



We provide evidence that the stripe instability results 
in a phase transition of first order, accompanied, in the 
thermodynamic limit, by jumps in various quantities, 
including the magnetization profile and the correlation 
length. This character of the transition seems to persist 
for vanishing pinning. 
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